Exact Zero Modes in Closed Systems of Interacting Fermions 
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We show that for closed finite sized systems with an odd number of real fermionic modes, even 
in the presence of many-body interactions, there are always at least two fermionic operators that 
commute with the Hamiltonian. There is a zero mode corresponding to the fermion parity operator, 
as shown by Akhmerov [l], as well as additional linearly independent zero modes, one of which 1) 
is continuously connected to the Majorana mode solution in the non-interacting limit, and 2) is less 
prone to decoherence when the system is opened to contact with an infinite bath. We also show that 
in the idealized situation where there are two or more well separated zero modes each associated 
with a finite number of interacting fermions at a localized vortex, these modes have non-Abelian 
Ising statistics under braiding. Furthermore the algebra of the zero mode operators makes them 
useful for fermionic quantum computation 
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Zero modes in non- interacting systems, i.e. eigenstates 
annihilated by a single-particle Hamiltonian, have a long his- 
tory in physics and in mathematics. Zero energy states are 
dissociated to certain types of topological defects in the back- 
ground fields in which electrons or quasiparticles propagate. 
The first example of such modes in physics appeared in the 
seminal work of Jackiw and Rebbi [3j] in one-dimensional 
and three-dimensional systems, where the topological defects 
were domain walls and hedgehogs, respectively. In both these 
Examples the physical consequence of the zero modes is the 
fractionalization of electron charge. Fractional charges can 
^Iso be bound to vortices in a Kekule dimerization pattern 
in two-dimensional graphene-like systems The zero mode 
solutions in two-dimensions were first found by Jackiw and 
Rossi in the study of Dirac fermions in the background 
of scalar and vector gauge fields of the Abelian Higgs model. 
In the condensed matter context this corresponds to a super- 
conductor (where charge cannot be fractionalized, since it is 
not conserved). The number of zero modes in such system 
of Dirac fermions in two-dimensions equals the magnitude of 
the net vorticity independent of the details of the profile of 
the Higgs fields, a result that was shown by Weinberg [6] to 
be tied to the index theorem. 

A modern example of a physical realization of the model 
m Ref.i was presented by Fu and Kane Q , who showed that 
a Dirac-type matrix equation governs surface excitations in a 
topological insulator in contact with an s-wave superconduc- 
tor. A vortex in the superconducting order parameter leads 
to a zero mode solution. Because of the reality conditions im- 
posed by the symmetries of the Bogoliubov-de Gennes (BdG) 
equations describing the superconductor within the mean- 
field approximation, the zero energy solutions correspond to 
Majorana zero modes, which are the focus of our study. Ma- 
jorana fermions are self-adjoint operators 7i which can be 
written as a sum of an annihilation and creation operator 
for one fermion mode and which satisfy the algebra: 



(1) 



[Hmf, li] ~ 0, these modes are in principle protected from 
decoherence as the mean field Hamiltonian, when restricted 
to the subspace generated by these modes, is zero. Recently 
it has been argued that quantum and classical fluctuations 
in open infinite systems (for example when the system is in 
contact to a bath) lead to decoherence of information stored 
in such modes Below, instead, we shall focus on closed, 
finite systems, which have markedly different properties from 
those coupled to an infinite environment. 

The purpose of this letter is to study zero modes of in- 
teracting many-body fermionic Hamiltonians, beyond mean- 
field approximations. We will assume that the relevant de- 
grees of freedom may be described by an odd number of 
Majorana fermions {71, 72, . . . , 72Af-i-i}- This formalism also 
handles the case when complex fermions are present, as 
we may change basis from complex to Majorana fermions: 
Cj = \ {j2j + «72j+i) ^c^j = k (72j - «72j+i)- For an interact- 
ing many-body Hamiltonian, a zero mode means a Hermitian 
fermionic operator 



(2) 



written as a multinomial with sums and products of 7i's, 
that commutes with the Hamiltonian, [H, O] = 0. For any 
such operator, O, exp (itH) Oexp (—itH) = O for all times 
t. As such there is no decoherence of the information stored 
in the correlators of such operators. 

We will find below, for systems of interacting fermions, 2^ 
linearly independent solutions of the form given in Eq. 
We will also extend our results to the case when interac- 
tions include bosonic modes (with finite dimensional Hilbert 
space) coupled to the Majorana modes. 

Quadratic Hamiltonians: Let us start, as a warm up, 

^J2^.n'^'^■J lil] with 



Gauss 



Because they are zero modes of some mean field Hamiltonian, 



with the simplest case where H 
hij =^ and hi^j real. We note that any quadratic Hamil- 
tonian may be written in this manner. Generic eigenoper- 
ator solutions satisfying , O \\ = \ 0\ are obtained 

by computing the commutators for operators of the form 
O — Q;i7i using the relations Eq. ([1]), and matching the 
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coefRcients multiplying each operator 7^ on both sides of the 
equation. One arrives in this manner at an eigenvalue equa- 
tion for the matrix 



/ /li,2 ^1,3 

h2,i ■■. 
h3,i ■■■ 



hi, 



2N+1 



\ 



2N+1,1 



l'2N+l,2N 



h2N,2N+l 





(3) 

The elements of the matrices T^Gauss ^^^^ ^ closely re- 
lated because the theory is Gaussian - there will be modifi- 
cations in the case of interacting systems. Note that T^^auss 
is an odd-dimensional Hermitian antisymmetric matrix so 
it has an eigenvector with zero eigenvalue and real com- 
ponents (ai, a2, ■ ■ ■ , a2N+i) which corresponds to the zero 
mode O = J^i ctili- Notice that it follows from the relations 
in Eq. (H]) that = O and = Ei x 1- 

Let us now introduce notation so as to arrive at the same 
■^Gauss ^ ^g^y t\ia± will be similar to the calculations 
for interacting systems below. Matching the coefficients 
multiplying each operator 7^ on both sides of the equation 
I^^Gauss^ _ ^ ]_|g achieved easily if we think of the 

7i as basis vectors and define an inner product for operators 
A and B as {A,B) = CoeS^{A^B), where 



CoeS^ z 1 + ^ 7i + ^ j^j 



(4) 



i.e., the function CoeS^{Q) returns the coefficient propor- 
tional to the identity in the multinomial expansion of the 
operator Q. One can check that the inner product is Hermi- 
tian, {A, B) — {B, A)* and it follows from the algebra of the 
7i's that the inner product gives (7^,7^) = Sij. 

Armed with this inner product we then compute the ma- 
trix 



Gauss 



(7.,[i^^^"^^7,]) 

-(7„[7?«^"-,7.]) 



-n 



Gauss 



(5) 



where the last line follows by direct computation and the 
fact that hij — —hj^i £ R. Once again Ti^^^^^ is given by 
Eq. ([3]) above. We thus arrive once more at the result that 
zero modes can be determined from null vectors of a linear 
eigenvector equation for a Hermitian anti-symmetric matrix 
^Gauss qJjJ dimension). 

Quartic Hamiltonian: We will consider a Hamiltonian 
given by: 



commute with TjQuart^ ^jjj work with a vector space 
that is spanned by all linearly independent Hermitian modes 
obtained from products of individual Majorana fermions 7^ : 



O7 
I7 
27's 
37's 

27V + 1 7's 



1, 

71:72,73, • • ■ ,72Af+l, 

ni72, i7l73, ■ • ■ , j72Af72Af-l-l, 

-1717273, . . . , — j72Af-l72Af72Af+l, 



(7) 



There are in total Y.TJ^ (''\+') 



■ 72Af+l • 

: 2^^+^ such operators, 
which we will denote by Ta, for a = 1, ... , 2^^+^. For each 
a we define tIq to be the number of 7's in the product Tq, 
and we let L{a) = {ii(a), . . . ,i„^(a)} be the list of indices 
appearing in the product Tq. With this notation, one can 
write 



,(na-l)/2 



7ii (0)7*2(0) • • • 7i„a (f ) ■ 



(8) 



The choice of phase factor guarantees that Tq ~ T\ and 
= 1. Using Eq. ([5]) one verifies that, up to a phase, the 
product of two Ta's gives a third: Ta = Tc, where 

c satisfies L(c) = L{a) U L{b) \ L{a) n L{b) and s (a, b) e N. 
Without loss of generality, we shall reserve the labels a = 1 
and a = 2^^+^ for the identity and the total parity operators: 
Ti = 1 and T22iv+i = i^'^^+^^'^jij2 ■ ■■I2N+1 = Tpcrmi- 
We can now rewrite the Hamiltonian Eq. © as 



a\n{a)—2 



ha Tn 



a\n(a)—A 



(9) 



for some coefficients ha , Va defined when 71(a) = 2 or 4, 
respectively, and ha, Va E K. Below we will convert i^Q^^rt 
into an operator acting on the vector space spanned by the 
Tq's with the action being given by the linear transformation 
where ijQ^art ^^^^g commutation: O [ijQuart^ ^^j ^ ^ 

first step we extend the inner product given in Eq. ^ above 
to the space spanned by T^ i.e. {A,B) = CoeS^{A^ B). 
One can check that the inner product is Hermitian, {A, B) = 
[B, A)* and the set Ta forms an orthonormal basis. Further- 
more, up to a multiplicative constant, we see that it is also 
given by the usual trace inner product: 

iA^B)^,^^tr{A^B). (10) 

Here, tr is taken over the space spanned by Ta. Indeed this 
can be checked by noting that Eq. (fTU)) is linear, so it is suf- 
ficient to consider only terms of the form A = T a, B = T^. 
There are two possibilities: 1) Ta = Tb in which case 
tr (Ta^Tfe) = 2^^+^ (the dimension of the vector space) 
2) Ta ^ Tb, for which case tr (Ta'''Tfc) = 0, and Eq. (fTOl) 



2J 



lilj ~^ j k I 'Ji^j'Jk'yi, (6) holds. Wc uow couiputc the matrix elements Since 

is an anti- Hermitian operator (or i times a Her- 
mitian operator) all the matrix elements of H^^"^'' are imag- 
inary. Now because {Tf,} is an orthonormal set we may com- 
pute matrix elements by taking inner products: 



with hij a real and anti-symmetric matrix and Vij,k.i real 
and antisymmetric under odd permutations of i,j,k,l (we 
have dropped an irrelevant constant that gives a state in- 
dependent energy shift). We will look for operators that 



n 



Quart 



(Ta,[i?Q""'-*,Tb]) 



(11) 
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so 'H^f,"'*'^* is antisymmetric. The equality in the last line of 
Eq. (fTT|) comes from the cyclic property of trace. Therefore 
we arrive at a Hermitian anti-symmetric matrix T^Q"'^"'*. So 
far, this matrix has dimension 2^^+^ x 2^^+^, which is even. 
However, one can break this matrix into four block-diagonal 
pieces. First, because TjQu^rt contains only even Tc, that is 
with Uc even, sectors with opposite parity are not mixed by 
^S)"'"^*' so necessarily Ua = nt, mod 2. Therefore we break 
^Quart jjj^Q blocks acting on the fermionic and bosonic {T^}, 
each block a 2^^ x 2^^ matrix. Second, notice that both the 
identity and the fermion parity operator commute trivially 
with ijQuait^ gQ they each reside in a 1 x 1 block. The identity 
is in the even sector (ni = 0) and the fermion parity operator 
is in the odd sector (npormi = 2N + 1). Therefore we have 
broken down T^Q"^''* into four odd-dimensional Hermitian 
and anti-symmetric block matrices: there are four operators 
that commute with the Hamiltonian //Quart _ ^j. ^gj-Q mode 
solutions. They are, in the even block, the trivial identity 
Ti = 1 and the Hamiltonian i/Q^art p].opg]._ ^nd in the odd 
sector the fermion parity Tpcrmi [1] and another non-trivial 
solution O — Oa Ta, with aa solutions of TL'^^^'^^ at = 
0. 

Generic Fermionic Hamiltonians: Let us al- 
low for arbitrarily high order interactions. That 
is we will consider Hamiltonians of the form 

JjGen ^ + T.^,jM^^^,3,k,ll^lJlkll + 

«E»j,fc,i,m,n<3»j,fc,;,™,«7«7j7fc7i7m7n + • ■ •, which may 
also be expressed as 



a\n{a) — 2 a|n(a)— 4 a\n{a)—6 

(12) 

where ha, Va, Qa,--- £ R- We can construct the matrix 
^Gcn gjmilarly to what we did above, it is still a Hermitian 
antisymmetric matrix. Nothing changes in the argument, 
and the essence is that the Hamiltonian contains only Tc 
with even Uc, and therefore one can break 'H'^''" into four 
block diagonal pieces exactly the same way we did for quartic 
Hamiltonians and obtain zero modes. 

Bosonic Modes: We now partially extend our ideas to the 
case of an odd number of Majorana fermions coupled to some 
bosonic modes. Our main limitation is that in order to in- 
sure convergence, to have finite dimensional matrices only - 
we will "truncate" the Hilbert space of the bosonic modes to 
a finite number of states. More precisely we will assume that 
the relevant Hilbert space for the bosons is M dimensional 
and labeled by the states {|1) , |2) ... \M)} As such we 
may represent all boson operators by M x Af Hermitian ma- 
trices. One can then write a Hamiltonian that generalizes 
Eq. dUI): 



H 



Gen— Bosc 



= e 



MxM 



a\n{a)—2 



MxM . 



System 



Environment 




Figure 1: The system in tunneling contact with the environment. 
The system is composed of CdGM states [llj, while the environ- 
ment is everything else. 



+ E ^a"^'"^ ® T„ + ^ Qf ^Ta + ... 

a\n{a)—4: a\n{a)—Q 
M^ 

= E E Wa,pTa<E>hp, (13) 

a\'n{a) even p— 1 

with e*^^*^ /if Va^'^""^, Qa""^^ Hermitian matrices 
and we expanded the bosonic M x M Hermitian matrices into 
an orthonormal basis {hi, h2, ...hjij-i}, with {hp,hq)Bosa = 
Spq. The inner product is {A, B)bosc = JT^^i^^^)- 
not too hard to see that this is a positive definite symmetric 
form on the space of bosonic operators [l3| • Without loss of 
generality, we take hi = ^mxm- 

We can combine the operators in the fermionic and bosonic 
spaces and define ^a,p = Ta(E)hp, with the usual tensor space 
inner product [l3|- These states are orthonormal because 

= S„,_h5r,_„. We can 



{^a,q, ^b,g) 



total 



(Ta,Tf,) X {hp,hg)i 



also check that this is expressible as a trace: {A, B)^^^.^^ = 
2iN+i jjtr (^^-S)- Here the trace is over the total space 
spanned by ^a,p- 

Armed with these combined operators, we can show that 
there is an exact zero mode in exactly the same way we have 
done in the previous case. We need the matrix: 



n 



Gen— Bosc 
a,p;b,q 



Gen — Bosc 



^a,p] ) 



-H 



(14) 

Gen— Bose 



which is Hermitian and anti-symmetric. The last equality 
in Eq. can be checked similarly to Eq. ([TT|l . We then 



break 'H^'j^i^f°^° into even and odd block diagonal spaces, as 
before. In this way, we find two zero modes in the even sector, 
Ti (g) ft-i = 1 (g) Imxm, and i/Gon-Boso proper, and two zero 
modes in the odd sector, Tpormi ® ImxM and another non- 
trivial solution O — J2a p'^O'^p'^'i ® with aa,p solutions 



= 0. 



Counting of modes and their structure: Let us count all 
zero modes in the system. We first start with the Gaussian 
part of the theory, including bosons, and then later we add 
the interactions. Consider a Hamiltonian given by: 



M 

^Gauss \ ^ 



m) TO + - 



1^ 



2 ^ 



«72i72j + l 



(15) 
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(Notice that i72j72j+i — '^c\ci — 1.) By inspection, there 
are M x 2^ bosonic zero modes ah given by operators 
of the form Olj,^^ ^ \m) (m| ® l\f=i (*72,72,+i)'^ with 
TO = 1, . . . ,M and 9j = 0, 1 for j = 1, . . . ,N. There are 
similarly M x 2^ fermionic zero modes, simply given by 
(i^Fcrmi ^ 0^°?° , 7, . These zero modes have a nice alge- 

braic structure: 1) they are all Hermitian, 2) all of them 



-jFcrmi/Bosc 



= 1, and 3) all zero modes 
= 0. As such any one 



square to one: {o'^-^ie,} 

commute- [(nF^'n'/Boso .^Fcrmi/ 

commuie. ^^„^^g.y ^'~-'m',{9'^ 
of the fermionic modes, and only one mode at a time, can be 
used for fermionic quantum computation 

Let us now show that the number of zero modes and 
their commutation relations do not change in the presence 
of weak interactions. To do so, as a first step, consider 
the following family of Hamiltonians H^^^ = i^Gauss ^ 
I^m,{e,} ^^foj} '^ith ^m,{e,} e K, and we note that 

{<^m,{ej}} G E^^^". It is not to hard to see that other 
then for points of accidental degeneracy all zero modes of 
all Hamiltonians of the form H^^^ are given by O^™'}^"*'"- 
As the next step, consider zero modes of Hamiltonians given 
by = f/t ijl-s} u. All the zero modes are now given 

by 0^™^^°^'^ U , and as such also satisfy conditions 1), 
2), and 3) of the previous paragraph. As before, exactly 
one mode from the fermionic set can be used for quantum 
computation Q. To complete the discussion of the count- 
ing and structure of the zero modes for interacting systems, 
it remains for us to show that any Hamiltonian with weak 
interactions can be written as a 

To show this, we consider the map T : U (A/^ x 2^^) © 
RMX2" ^ rM^x2^" gi^gj^ hyT{U, {S,n,{e,}}) = WH^^^U. 
It is enough to show that the image of U (^AP x 2^^) © 
■[^Mx2 contains a small open neighborhood of TjGauss^ in- 
deed, as any sufficiently weakly interacting Hamiltonian 
can be found in a small neighborhood of a non-interacting 
one this would show that U^H^^^U is a representation 
of all sufficiently weakly interacting Hamiltonians. By 
the implicit function theorem it is enough to show that 
dT is a surjective mapping onto R*^ . Now writing 



U = 



we get dT(H,{S^,{g^y}) = 



^ ^Gauss 



+ 



J2m,{ej} ^rn,{ej} 0^,fe.}- From this we see that all the zero 
modes are explicitly in the image of dJ-. Since the transfor- 
mation * — i *, is an invertible linear operator 

when restricted to the space of all non-zero modes, all non- 
zero modes are also in the image of dT as well. As such all of 
gAf x2 -j-j^g image of dJ-. This shows that up to conju- 

gation by a unitary transformation the structure of the zero 
modes is the same as in the non-interacting case completing 
the proof. 

Comparison with previous work: In Ref. [1], the fermion 
parity operator Tpormi was discussed. This parity operator 
commutes with any Hamiltonian, since it is formed by the 
product of all the operators ji. This operator sits on its own 
1x1 block of the matrix H, for all cases studied, including 



in our generalization that includes bosons interacting with 
the fermionic modes. 

In contrast, the other zero mode solutions found in the 
larger odd-dimensional block of H do depend on the form of 
the Hamiltonian. There are M x 2^ — 1 of them. Furthermore 
one of the modes has a particularly simple structure O = 
e^^ X^i '^ili e~'^ which is continuously connected to the non 
interacting mode (consider Ot = e**^ 0^7^ e~**^). This 
mode is different from the fermion parity mode [l| and, as we 
shall see below, for weak interactions (small H) it is better 
protected from various forms of decoherence when the system 
is coupled to a generic bath. 

Decoherence: Consider the setup shown in Fig. ([T]). We 
consider a simple perturbing tunneling Hamiltonian of the 
form: AH = i U^fiiji, with ti G M. Here rji refer to Hermi- 
tian fermionic modes relevant to the environment. In previ- 
ous works it was demonstrated that (0 (0) O (T)) is a good 
measure of the coherence of a qubit composed of localized 
Majorana modes Here O is an operator used to encode 
the qubit, and we will assume that the qubit and environ- 
ment start uncorrelated. By Taylor expanding e^"^^^ and 
keeping only leading order terms we obtain {O (0) O [T]) = 

1 - \T^Y.,,Uh {M X {(07.7,0) + (O7.O7,)} 

+ (77,77.) x { (07,07.) + (OH-70}}- 

(16) 

We can understand how this expression scales for vari- 
ous operators, in particular for O = Tq, tIq odd, we get 
that (T,(0)T,(r)) = l-2T^Y.^eL(a)tHv^)^^^- Since 

('^^Env — operators with larger Ua decohere more 
quickly, at least for short times. This indicates enhanced 
stability for operators that are similar to single Majorana 
fermions, like the new zero modes presented here. 

Braiding: We would like to consider the idealized 
case of several Fermi zero modes {O^}, of the form 
e*^* Q^i7i e~*^* , each corresponding to its own individual 
finite environment and labeled by t. We further assume that 
the individual environments do not interact with the rest of 
the system. In this case because the modes are composed of 
sums of products of odd numbers of Majorana modes we see, 
following Ivanov [12], that the transformation properties of 
two zero modes under exchange are given by: 

01 -> O2 (17) 

02 ^ -Oi. 

The minus sign comes when vortex 1 crosses the "cut" cor- 
responding to vortex 2. These rules are identical to Ising 
braiding rules. To extend the derivation of Eq. ([T7| given 
in |12], we must show that the many body holomony when 
restricted to the zero modes is zero. Indeed referring to 
^13ril5i] , we know that in order to add in the many body 
holomony, in the Schrodinger picture, it is sufficient to con- 
sider Hamiltonian evolution within the zero energy subspace; 
with a Hamiltonian whose matrix elements are given by 
-fff2,o' — *(^l^l^')- Here jfi) and \VL') are instantaneous 
zero energy eigenkets. In the Heisenberg picture, this evolu- 
tion corresponds to an evolution of the operators Ta when 



5 



acted on by the Hamiltonian 'Ha,p-b,q, see the discussion fol- 
lowing Eq. ([9]). As such it is enough to show that any Hamil- 
tonian 'Ha.p;b,q cau Only have zero matrix elements within 
the subspace of zero modes. Indeed we know that up to a 
unitary transformation the structure of the zero modes is 
the same as in the non-interacting case, so by transform- 
ing H U^HU = H we reduce the problem to the non- 
interacting case. As such it is enough to show that any 

Ha.p-b.q — when restricted to the space of zero modes. 
By linearity it is enough to consider only Hamiltonians of 
the form \m) {n\ Tfc. By explicitly taking commutators 
with Om.{0j} we see that all matrix elements within the 

zero energy subspace are zero {Ha,p-b,q = 0). This means 
that under braiding the fermionic zero modes transforms as 

Ci,m,{ej} C'2,m,{ej}, Oi,m,{ej} -C2.m,{0,} fis'l • In par- 
ticular the non-interacting holomony, Eq. (|17ll . is recovered. 

Conclusions: We presented a systematic treatment of 
closed interacting systems with an odd number of real 
fermions. This formulation allowed us to find the zero mode 
solutions of interacting Hamiltonians, i.e., operators that 
commute with the many-body Hamiltonian. In addition to 
the fermion parity operator that can be viewed as a con- 
stant of the motion for any Hamiltonian, we have found the 
solution that connects continuously to the Majorana mode 
for non-interacting systems as the interactions are switched 
off. These modes couple more weakly than the fermion par- 
ity mode to an environment once the system is opened up 
to an outside infinite bath [8]. Therefore, the solutions that 
are continuously connected to the non-interacting Majorana 



modes should lead to slower decay rates in the presence of 
a bath. We have also verified that, under idealized condi- 
tions when multiple such modes exist, they obey Ising like 
statistics under braiding. 
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